Electron properties of Carbon nanotubes in the field effect regime 
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Electron properties of Carbon nanotubes can change qualitatively in a transverse electric field. 
In metallic tubes the sign of Fermi velocity can be reversed in a sufficiently strong field, while in 
semiconducting tubes the effective mass can change sign. These changes in the spectrum manifest 
themselves in a breakup of the Fermi surface and in the energy gap suppression, respectively. The 
effect is controlled by the field inside the tube which is screened due to the polarization induced on 
the tube. The theory of screening links it with the chiral anomaly for ID fermions and obtains a 
universal screening function determined solely by the Carbon tt electron band. 



The possibility to change electron spectrum of Carbon 
nanotubes by external field is of interest for basic research 
as well as for nanoscale device engineering. Carbon nan- 
otube (NT) is a \D metal or semiconductor depending on 
the chiral angle Metallic behavior can be suppressed 
by parallel magnetic field that induces a minigap at the 
band crossing Similar minigaps appear in nominally 
metallic chiral NTs due to the intrinsic curvature 
Novel properties are predicted for the BN tubes having 
no inversion symmetry . 

Here we examine the changes induced in the NT elec- 
tron spectrum by transverse electric field £ strong enough 
to mix different NT subbands: 

e£R~ ^n = hv/R, f [MV/cm] ~5.26/i?2[nm2], (1) 

where R is the tube radius and v is electron velocity. In 
such a field the effect on electron spectrum is dramatic: 
in metallic tubes the electron velocity v — de/dp can be 
reduced and even reverse the sign, causing Fermi sur- 
face breakup, while in semiconducting tubes the effective 
mass sign can change, which is accompanied by strong 
suppression of the excitation gap (Fig.|l]). 

The NT electron system in this regime can be a host of 
intriguing many-body phenomena. The reduction of elec- 
tron velocity in metallic tubes leads to an increase of the 
dimensionless interaction e'^/hv that controls the Lut- 
tinger liquid properties [^. One expects this to enhance 
the intrinsic, interaction-induced energy gap predicted to 
be small in pristine NTs |^ . 

Even more peculiar is the negative v state with in- 
tertwining electron and hole Fermi surfaces (Fig.|l| top). 
This system provides a realization of a metallic state un- 
stable with respect to electron-hole pairing into excitons. 
Such an instability, long-envisioned |^ by Mott, Kcldysh 
and Kopaev, and others, is especially interesting for the 
mirror symmetric electron and hole bands, described by 
an analog of the BCS theory. The chiral gauge symmetry 
of the NT electron Hamiltonian discussed below elim- 
inates the interband matrix elements of particle density. 
This makes the phase of the excitonic order parameter a 
gapless Goldstone field, similar to the BCS order parame- 
ter. Despite its simplicity, no system with such properties 
has been unambiguously identified so far, and it is thus 



possible that nanotubes in a transverse field present a 
unique opportunity to study this phenomenon. 

The field (|l|) required to create this state is to be 
achieved within the tube where the external field is 
partially screened However, because of the dis- 

crete bands with relatively large interband separation 
Aq [eV] = hv/R = 0.53/i?[nm], the transverse field pen- 
etrates in the NT fairly well. We find the screening factor 
to be i?- independent and close to 5 (in accord with |^ ) , 
for both metallic and semiconducting tubes. With the 
screening taken into account the numbers for the required 
field remain feasible. (Fields up to 50 MV/cm have been 
demonstrated recently in 2D gated structures ||ll|] .) 

Unexpectedly, the problem of transverse field screening 
has a relation with the chiral anomaly. A significant part 
of electron energy in an external field, in the massless 
Dirac model (standard for NT), arises from the effects 
at the Fermi sea bottom, where a regularization of this 
model is required. However, the anomaly links the reg- 
ularized energy with the properties near the Fermi level 
(the number of fermion species and their velocity), and 
thereby generates a universal screening function deter- 
mined only by the Carbon tt electron band. 

We consider the NT at half-filling ignoring charge accu- 
mulation due to gating. Gating in itself will not modify 
the transverse field within the tube, since a uniformly 
charged cylinder is equipotential. Charging may affect 
the inner NT field indirectly via changing screening, but 
this effect should not be significant at moderate gating. 

The TT electron Carbon band is described by a nearest- 
neighbor tight-binding problem eipr = t ipr' with the 
hopping amplitude t ~ 3eV. The states with small 
e| <C i are described, separately at each of the K and 
K' points, by a massless Dirac Hamiltonian 



V = 3tac^c/2h (2) 



For a NT in the presence of a transverse electric field £, 

Ho — hv {iaxdy + ayk) — e£ R cos{y / R) (3) 

with k the longitudinal momentum and R the NT radius. 
The boundary conditions are quasiperiodic: 



, 



metallic 
semic. 



(4) 
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The effects of NT curvature g as well as of a parallel 
magnetic field can in be included by slightly shifting 
S away from the ideal values . 

We employ a chiral gauge transformation 



e£R^ 

hv 



My/R) (5) 



which preserves the condition and turns Eq.(||) into 

Ho = hv (^ia^dy + fce2*'^=='^(^Vj,) (6) 



The transformed Hamiltonian reveals that, in particu- 
lar, the spectrum at fc = is totally independent of the 
transverse field (FigM. 
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FIG. 1. Electron bands transformation: velocity reversal 
in metallic NT (top); effective mass sign change in semicon- 
ducting NT (bottom). The bands are shown for the dimen- 
sionless transverse field u = e£R/Ao below and above critical. 
Inset: Energy gap suppression in a semiconducting NT. 

The spectrum of a metallic NT near the band crossing 
at kR ^ 1 can be found by employing the degenerate 



perturbation theory. We project the Hamiltonian (|^) on 
the two states | T)i I i) degenerate at fc = 0, and obtain 

e(fc) = ± [nvjQ{2u)] k , u = eSR/Ao = eSR^/Tiv (7) 

Electron velocity changes sign at the roots of the Bessel 
function Jq, first at 2m = /ii « 2.405 (Fig.|l|). At u above 
critical the Fermi surface fractures: an additional small 
pocket appears for each spectral branch. 

The level shifts in semiconducting NT at small k are 
given by the second order perturbation theory in the k- 
term of the transformed Hamiltonian (|^): 



= (fc) -±Ao((n + 5) + A„(fci?)2) 



2(71 + S) 



4(n -|- (5)2 — 



(8) 
(9) 



For 5 = 1/3 the curvature of the lowest band Aq changes 
sign at Uc « 0.6215 (Fig.|]). This leads to a singular be- 
havior of the excitation gap which is constant at u < Uc 
and sharply decreases at u > Uc (Fig.|l| inset). This oc- 
curs because of the lowest excitation energy shifting from 
fc = at u < Uc to fc ^ at M > lie. The threshold- like 
suppression of the gap can be detected by a transport 
measurement in a thermally activated regime. 

The chiral gauge symmetry (||) that protects the spec- 
trum at fc = is a distinct feature of the Dirac model 
(^, (^. The TT electron tight-binding problem, in the 
next-lowest gradient order, generates a correction to the 
Hamiltonian (||) violating the symmetry (|^): 



H = no + Ae-^«"^[a, (dl-dl) -2ayd,dy] 



3ii 
e 2 



(10) 



with X = 



hv and e the NT chiral angle IJ. The 
transformation (||) applied to ( p^ ) gives a minigap |Q 



|sin36l| (a^^jShv) [eERf 



ill) 



Since 8i?, 3> Oc-c, the minigap ( |ll|) is too small to alter 
the behavior at the energies of interest, e ~ Aq. 

The main effect of electron interaction is screening of 
the inner field that couples to the electron motion. Here 
we derive the relation between the inner and outer fields. 
We first show how the screening problem is reduced to 
the calculation of the NT electron energy in the presence 
of an external field. Hereafter we measure all energies 
in the units of Aq = hv/R and use dimensionless field 
u = e£R/Ao. From Gauss' law, the fields inside and out- 
side the tube are related with the induced surface charge 
density per one fermion species (spin and valley) by 



£ 



An ■ 



(12) 



where the factor 1/2 accounts for depolarization in the 
cylindrical geometry. In Eq.(|l^ we projected the actual 
charge density on the cos ip harmonic as a{ip) —>■ Aa cos ip, 
ignoring the higher order harmonics. (Here ip = y/R.) 

To obtain the cos ip harmonic of the induced charge, 
we evaluate the dipole moment per unit length as P = 



2 



—dW{£)/d£, where W{£) is the energy of one fermion 
species as a function of the inner field. Combining this 
with the relation a — P/{'kB?) and with the Gauss' law 
(|l^), and passing to dimcnsionless Uoxt, u, we obtain 



Mext 



w + 8— P(u) 
nv 



(13) 



After the dipole moment P{u) is known Eq.(|13|) can be 
solved for the inner field u in terms of the outer field Uoxt- 
We consider the general problem of electron energy in 
a transverse field in a free particle model. The electron 
levels e„_fc perturbed by the field can be easily found nu- 
merically at each value of the longitudinal momentum k 
by using a transfer matrix for Eq.(^. The level shifts 

Se-n.k — en,fc(u) — \ decrease at large |n|, and the series 
the total change of the occupied states energy 

+00 

Eo{k)= Sen,k (such that e„,fc(u) < 0) (14) 

n— — 00 

rapidly converge at n — > ±(X). 

There are two basic problems with Eq.(p^: 1) Due to 
an upward shift of the filled levels (Fig.l), Eq is positive 
and also has positive derivative dEo/du. Hence Eq.(|l^) 
leads to the dipole moment Pq — —dEo/du opposite to 
the field, i.e. to an unphysical "diamagnetic" polariza- 
tion sign instead of the expected "paramagnetic" effect. 
2) The dependence of the energy Eq on the longitudinal 
wavevector k leads to an ultraviolet divergence in the in- 
tegral P ~ J P{k)dk, because Eo{k) increases with |fc|, 
saturating at |fc|i? 3> 1 at an asymptotic value ^v? ■ 

Both difficulties are resolved by taking into account a 
fundamentally important contribution to the energy that 
arises due to the effects at the Fermi sea bottom. Physi- 
cally, the finite electron band width invalidates the mass- 
less Dirac approximation at large negative energies. This 
contribution, however, depends solely on the number of 
Dirac fermion species and their velocity v, and is totally 
insensitive to any other details including the longitudinal 
momentum k value. We find that 



E, 



(15) 



for each fermion species. Remarkably, Eq.(|T^) can be ob- 
tained without detailed discussion of the behavior at the 
interatomic length scales — the universality of Eq. (|l5|) is 
rooted in the physics of the chiral anomaly in the 14-1 
fermion problem. The resulting total energy integral 



W 



(Eoik) + E^nom) 



dk 
2ir 



(16) 



converges after Eo{k) is offset by E'anom (Fig.||). 

We evaluate the energy (|l^ ) and derive the anomaly 
(|l|) for a weak field u < 1. The NT bands at u = are 



In a half-filled system with just the e„ (fc)<0 bands filled, 
the external field u changes the Fermi sea energy by 



dk 
2-K 



(18) 



Here the superscript in indicates regularization by 
truncating the interaction with the external field at a 
certain large negative energy. We check that this con- 
tribution to the energy is independent of the details of 
truncation and obtain the anomaly (^5|) by choosing a 
convenient truncation scheme. 
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FIG. 2. Dipole moment P(fc) = -d{Eo{k) + E^noi-n)/du 
per one fermion species in a semiconducting NT as a function 
of k. Note that the energy anomaly (^) cancels with Eo{k) 
at kR S> 1, assuring convergence of Ptotai ~ J P{k)dk/2TT. 
Note also that P{k — > 0) is dominated by the anomaly, since 
-Bo = at fc = due to the chiral gauge invariance 

The level shifts Se~{k), in the second order of the per- 
turbation theory in the external field V=—e£Rcos(p, are 



5e„ 



E 



\{m+\V\n-)\^ 



E 



|(TO-|y|?l-)|^ 



(19) 



;(fc) = ±v/(n + 5)2 + fc2 



— 00 < 71 < 



(17) 



where the superscript ± indicates the electron and hole 
branches and the k dependence is suppressed. Due to the 
integration over ip with V oc cos tp in the matrix elements 
the only nonzero terms in (^|) are those with m = n ± 1 . 

We now show that the sums over and in (p^), 
respectively, give the regular and the anomalous contri- 
butions to the total energy W — J2'^ Se:^. Different be- 
havior of the two sums under regularization stems from 
their different convergence type. Individual terms in the 
sum over e+ decrease rapidly at large m, so that the se- 
ries for W is absolutely convergent. On the other hand, 
in the sum over e~ the terms do not change at large m 
and thus the corresponding contribution to W is given 
by poorly convergent and regularization-sensitive series. 

Taking from (19) just the terms with e^j, evaluating the 
matrix elements (m+|y|n^) and summing over n yields 



3 



Eo{k) = 



n 



,{k)-{n + S)in'+S)-e 



(20) 



with n' = n + 1. The sum (^0|) rapidly converges at large 
n — > ±00 and can be easily evaluated numerically. 

Now we consider the sum of the level shifts W — 



'^'n^^n taking into account only the second term in (19). 
At the first sight this sum is identically zero. Indeed, 
due to the symmetry {m~\V\n^) — {n^\V\m~), in the 
sum over n with m = n ± I all the terms cancel in pairs. 
However, truncation of the interaction at a large negative 
energy compromises the cancellation and yields a finite 
result. If one sets {m~\V\n~) = for all |m| or |n| ex- 
ceeding a large number N, there will be just two terms 
in the sum over n that do not cancel: 



\{N'-\V\N-)\^ \{-N'-\V\^N-)\'' 



(21) 



with N' = N — 1. Evaluating the matrix elements and 
energy levels is straightforward because at large TV ^ |fc| 
one can set k — 0. The result, coinciding with (^5|), is 
robust under a change of the regularization. 

The expression (|l5| ) for the energy anomaly, derived 
above for the weak field, is in fact more general. To illus- 
trate this we consider a special case of zero longitudinal 
momentum fc = and derive (|l5|) from bosonization, 
without using perturbation theory in u ^ 1. After the 
problem (0) is bosonized in the standard way [12| , using 



we obtain a quadratic Hamiltonian 



n 



E [tid.^jf + td^^jU{y)\dy (22) 

"'0 , r o 



j=L,R 



The second term in ( p2|) representing interaction with 
the external field U{y) can be decoupled by a shift 
4>j cj)j — j-^ U{y')dy'. The Hamiltonian for (jy'j takes 
the form (02) with U — 0, while the ground state energy 



SE = 



2lTflVF 



2-rrR 



U'{y)dy 



(23) 



is nothing but the anomaly ( [Tsl) scaled by Aq — Tiv/R. 

After adding the energies ( po[) and (|l^), and integrat- 
ing in ([l6|) over k numerically, one obtains 



W 



-—u 
2 



0.196.. 
0.179.. 



for 5 = 1/3 
for (5 = 



(24) 



Eq.(|l2|) with the dipole moment P = —dW/dS and the 
charge density a = P/ttR^ yield the screening function 



hv ' 



(25) 



With e^/hv = 2.7 this gives S^^t/S = 5.24 for (5 = 1/3, 
and £cxt/£ = 4.87 for S — 0. The outer-to-inner field 
ratio ~ 5 (see [[lO[ for another derivation) renders the re- 
quired fields (|l|) feasible |jl^. Interestingly, the screening 
( 2^) is independent of the tube radius R and is almost the 



same in the metallic and semiconducting NTs. The latter 
is not surprising, since the screening is absent in a sin- 
gle ID mode approximation: the polarizability is related 
with dipolar transitions between different subbands. 

The radius-independence of (|2^) resembles an effect of 
a dielectric constant. We note, however, that the change 
of the inner field due to individual Carbon atoms polariz- 
ability is small in ac_c/27ri? -C 1. The result ( p5|) refiects 
the semimetallic character of the tt electron band with 
the density of states vanishing at the band center. 

In summary, nanotube electron states undergo inter- 
esting transformations in the field effect regime, leading 
to novel phenomena in both the single particle and many- 
body properties. The analysis of screening, performed 
using a relation with the theory of chiral anomaly, indi- 
cates that the fields required for the observation of the 
proposed effects are in the experimentally feasible range. 
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